New universal conductance fluctuation of mesoscopic systems in the crossover regime 

from metal to insulator 
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We report a theoretical investigation on conductance fluctuation of mesoscopic systems. Exten- 
sive numerical simulations on quasi-one dimensional, two dimensional, and quantum dot systems 
with different symmetries (COE, CUE, and CSE) indicate that the conductance fluctuation can 
reach a new universal value in the crossover regime for systems with CUE and CSE symmetries. 
The conductance fluctuation and higher order moments vs average conductance were found to be 
universal functions from diffusive to localized regimes that depend only on the dimensionality and 
symmetry. The numerical solution of DMPK equation agrees with our result in quasi-one dimension. 
Our numerical results in two dimensions suggest that this new universal conductance fluctuation is 
related to the metal-insulator transition. 
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One of the most important features in mesoscopic sys- 
tem is that the conductance in the diffusive regime ex- 
hibits universal features with a universal conductance 
fluctuation (UCF) that depends only on the dimension- 
ality and the symmetry of the system^ There exist three 
ensembles or symmetries according to the random ma- 
trix theory (RMT)2: (1). Circular orthogonal ensem- 
ble (COE) (characterized by symmetry index (3 = 1) 
when the time-reversal and spin-rotation symmetries are 
present. (2). Circular unitary ensemble (CUE) ((3 = 2) if 
time-reversal symmetry is broken. (3). Circular symplcc- 
tic ensemble (CSE) ((3 = 4) if the spin-rotation symme- 
try is broken while time-reversal symmetry is maintained. 
In the diffusive regime, the UCF is given by rms(G) = 
Cd/y/j3e 2 /h where a = 0.73,0.86,0.70 for quasi-one di- 
mension (ID), two dimensions (2D) and quantum dot 
(QD) systems and (3 = 1,2,4^ Although the RMT 
can apply to both diffusive and localized regimes, so far 
the universal conductance fluctuation has been addressed 
and established only in the diffusive regime. When the 
system is away from the diffusive regime, some universal 
behaviors have been observed. For instance, the conduc- 
tance distribution of quasi- ID systems ((3 = 1) with sur- 
face roughness was found to be universal in the crossover 
regime, independent of details of the system^ For quasi- 
1D systems with (3 = 1 , 2 the conductance distribution 
obtained from tight-binding model agrees with the nu- 
merical solution of DMPK equation^ In the localized 
regime, the conductance distribution of the quasi-lD sys- 
tem obeys log normal distribution. 5 In high dimensions, 
conductance distribution at the mobility edge of metal- 
insulator transition was also shown numerically to be uni- 
versal for 2D systems with (3 = 2, 4 and a 3D system 
with (3 = 1& In the localized regime, the conductance 
distribution of 3D systems is qualitatively different from 
that of quasi-lD systems^ It would be interesting to fur- 
ther explore the universal behaviors of these systems and 
ask following questions: Is there any universal behaviors 
away from mobility edge? Is it possible to have a UCF 
beyond the diffusive regime? If there is, what is the na- 



ture of the new UCF? It is the purpose of this work to 
investigate these issues. 

To do this, we have carried out extensive numeri- 
cal calculations for conductance fluctuations rms(G) in 
quasi-lD, 2D and QD systems for different symmetries: 
(3 = 1,2, 4. Our results can be summarized as follows. 

(1) . From diffusive to localized regimes, the conduc- 
tance fluctuation and higher order moments were found 
to be universal functions of the average conductance 
for quasi-lD, 2D, and QD systems and for (3 = 1,2,4. 

(2) . We found that there exists a second UCF near 
the localized regime for (3 = 2,4 but not for (3 = 1. 
Our results show that the new UCF depends weakly on 
the symmetries of the system and assumes the follow- 
ing value: rms(G) = Cd,e 2 /h. Here for (3 = 2 we found 
c d = 0.56 ± 0.01, 0.68 ± 0.01, 0.58 ± 0.01 for quasi-lD, 2D 
and QD systems, respectively while for (3 = 4 we have 
d d = 0.55±0.01,0.66±0.02,0.56±0.01. The conductance 
distribution in this new regime was found to be one-sided 
log-normal in agreement with previous results^ (3). In 
the localized regime with (G) < 0.3, the conductance dis- 
tribution does not seem to depend on dimensionality and 
symmetry. (4). For quasi-lD systems, the numerical so- 
lution of DMPK equation- agrees with our results. (5). 
For quasi-lD systems, the new UCF occurs when the lo- 
calization length £ is approximately equal to the system 
size L, i.e., £ ~ L for (3 = 2, 4. For 2D systems, we found 
that the new UCF occurs in the vicinity of the critical 
region of metal-insulator transition. 

In the numerical calculations, we used the same tight- 
binding Hamiltonian as that of Ref[^. Static Anderson- 
type disorder is added to the on-site energy with a uni- 
form distribution in the interval [—W/2,W/2] where W 
characterizes the strength of the disorder. The conduc- 
tance G is calculated from the Landauer-Buttiker for- 
mula G = (2e 2 /h)T with T the transmission coefficient. 
The conductance fluctuation is defined as rms(G) = 



(G 2 ) — (G) , where (• • •) denotes averaging over an en- 
semble of samples with different disorder configurations 
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FIG. 1: (Color online) Average conductance (a,b,c) and its 
fluctuation (d,e,f) vs disorder strength W for different sym- 
metry index (5 in quasi- ID systems. Insets: localization length 
vs W for p = 2, 4. 



of the same strength W. In the following the average 
conductance and its fluctuation are measured in unit of 
e 2 /h; the magnetic field is measured using magnetic flux 
4> = fisB/t where /is is Bohr magneton and t is the 
hopping energy that is used as the unit of energy. 

We first examine average conductances (G) and their 
fluctuations rms(G) vs disorder strength W in quasi- 
ID systems with different symmetry index (3 (see Fig. I). 
In our numerical simulation, the size of quasi-lD sys- 
tems are chosen to be 40 x 2000 for /3 = 1,2 (Fig. la, b) 
and 40 x 800 for (3 — 4 (Fig.lc). Each point in Fig.l 
is obtained by averaging 9000 configurations for /3 = 1 
and 15000 configurations for (3 = 2,4. In Fig.l, data 
with different parameters are shown. For instance with 
(3 = 2, (G) and rms(G) vs W are plotted for differ- 
ent Fermi energies with fixed <f> = 0.0628 and differ- 
ent magnetic flux with fixed Fermi energy E = 3.8. 
From Fig.l, we see that in the diffusive regime where 
(G) > 1 there is a plateau region for rms(G) with the 
plateau value approximately equal to the known UCF 
values rms{G) = 0.73/^3 = 0.73,0.52,0.365 (marked 
by solid lines). This suggests that one way to identify 
UCF is to locate the plateau region in the plot of rms{G) 
vs disorder strength and the plateau value should cor- 
respond to UCF. This method has been used to iden- 
tify universal spin-Hall conductance fluctuation^ that 
was later confirmed by RMT. 11 Importantly, there ex- 
ists a second plateau region for (3 — 2,4 but not for 
(3 = 1. The new plateau value approximately equals to 
rms{G) = 0.56±0.01 for f3 = 2 and rms(G) = 0.55±0.01 
for f3 = 4. In this regime, we found that (G) < 1 which 
clearly corresponds to the crossover regime. To confirm 
that the first and second plateaus are indeed in the dif- 
fusive and crossover regimes, respectively, we have cal- 
culated the localization length £ of the quasi-lD system 



FIG. 2: (Color online) Conductance fluctuation vs average 
conductance for quasi- ID systems (a,b), QD systems (c), and 
2D systems (d). 



(insets of Fig.l). It is clear that near the first plateau 
where W ~ 0.4 for (3 = 2 and W ~ 1 for (3 = 4 we have 
£ 3> L with L the length of quasi-lD system while near 
the second plateau we have £ ~ L (see insets of Fig.l). 

According to the UCF in the diffusive regime, it is 
tempting to conclude that this second plateau should cor- 
respond to a new UCF. However, in making such a claim 
one has to answer following questions: (1). whether the 
plateau behavior is universal? (2) . if it is, whether it can 
be observed in a wide range of parameters? (3). how 
is our result compared with the theoretical predictions 
whenever available? (4). whether such a universal behav- 
ior exists in high dimensions? In the following, we pro- 
vide evidences that the new plateau indeed corresponds 
to a new UCF. 

To answer the first question, we plot rms(G) vs (G) 
in Fig. 2a by eliminating W. The fact that all curves 
shown in Fig.l with different parameters (Fermi energy 
E, magnetic flux <J), and SOI strength t so ) collapse into 
a single curve for each (3 strongly indicates that rms(G) 
vs (G) is universal. To further demonstrate this univer- 
sal behavior, we have calculated the conductance fluc- 
tuation for a quasi- ID system with both magnetic flux 
and Rashba SOI. Although the Hamiltonian of this sys- 
tem is still unitary, both time reversal and spin rotation 
symmetries are broken. According to the diagrammatic 
perturbation theory^, the UCF is reduced by a factor 
of 2 when SOI is turned on. From RMT point of view, 
both systems (B ^ 0, t so = 0) and (B ^ 0, t so £ 0) 
are unitary ensembles and obey the same statistics. The 
fact that energy spectrum for (B ^ 0, t so = 0) is dou- 
bly degenerate accounts for the factor of 2 reduction for 
the system (B ^ 0, t so ^ 0). In Fig. 2a, we have plotted 
rms(G) vs (G) for the system with (B ^ 0, t so ^ 0). 
Once again, we see that all data from different parame- 
ters collapse into a single curve. If we multiply this curve 
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FIG. 3: (Color online) Average conductance and its fluctu- 
ation vs s for quasi- ID systems [(a) and (c)] and compared 
with results of DMPK equation for different N — 4,6, 18 
[(b) and (d)]. The arrow points the direction of increasing TV. 



by a factor of 2, it collapses with the curve of P = 2 (see 
inset of Fig. 2a). 

From Fig. 2a we see that in the diffusive regime with 
(G) > 1 the long plateau for each (3 corresponds to the 
known UCF (marked by solid lines). For P = 1 there 
is only one plateau. For (3 — 2,4, however, a second 
plateau region is in the neighborhood of (G) ~ 1. For 
[3 = 2, rms{G) is approximately a constant hence uni- 
versal in the region (G) = (0.6, 1.4) while for (3 = 4 this 
region narrower with (G) = (0.7, 0.9). Looking at Fig. 2a, 
it seems that the second plateau region is narrower than 
the first one. But if we look at rms(G) vs W where W 
can be controlled experimentally, the crossover region is 
enlarged since in the crossover regime (G) is not very 
sensitive to W while in the diffusive regime it is the op- 
posite. Indeed, in Fig. If, we do see that the ranges of 
the first and second plateau regions are comparable. If 
we fix W and plot (G) and rms(G) vs the length of the 
system, the crossover regime is enlarged further. These 
results are shown in Fig. 3a, c where the symbols represent 
our numerical result and solid lines correspond to exact 
solution of DMPK equation (to be discussed below) . We 
see that the window of the second plateau is much larger 
than the first one. 

Since the statistics of transmission eigen-channels of 
quasi-lD systems can be described by DMPK equation, 
we have numerically solved it- 3 for (3 — 2, A and com- 
pared our numerical results of tight-binding model with 
that of DMPK. Fig. 3b and Fig. 3d show the numerical 
solution of DMPK equation for (3 = 2, 4 with different 
N = 4, 6, ...,18 where N is the number of transmis- 
sion channels. Here s — L/l where L is the length of 
quasi-lD systems and I is the average mean free pathi^^ 
Fig. 3 clearly shows that in the diffusive regime where 
1 <C s <C N the first plateau corresponds to the usual 



FIG. 4: (Color online) Average conductance and its fluctua- 
tion vs W for QD systems with /3 = 1, 2, 4. 



UCF and there exists a much wider second plateau in 
the crossover regime where s and N are comparable with 
plateau value equal to our newly identified UCF. Fig. 3a 
and Fig. 3c show the comparison between the results of 
DMPK and that of quasi-lD tight binding model. The 
rms(G) vs (G) of DMPK equation is plotted in Fig.2b 
where selected data from Fig. 2a is also plotted for com- 
parison. The agreement between numerical and theoret- 
ical results is clearly seen. 

Now we examine the conductance fluctuation for QD 
systems. In the numerical calculation, the size of QD is 
L x L with L = 100 and the width of the lead L = 10 
for j3 = 2,4 while for (3 = 1 we used L = 150. Fig.4 
depicts (G) and rms(G) vs W for p = 1, 2, 4. Each point 
in Fig.4 was obtained by averaging 9000 configurations 
for P = 1 and 20000 configurations for /3 = 2,4. Similar 
to quasi- ID systems, we see only one plateau region in 
the diffusive regime for P = 1 with plateau value close 
to the known UCF value rms(G) — 0.70. In addition 
to the first plateau in the diffusive regime, there exists 
a second plateau for P = 2,4 which we identify to be 
the regime for new UCF. The new UCF is again in the 
crossover regime where (G) ~ 1 with the value UCF(/3 = 
2) = 0.58 ± 0.01 and UCF(/3 = 4) = 0.56 ± 0.01. In 
Fig. 2c we plot rms(G) vs (G). It shows that all curves 
for each P collapse into a single curve showing universal 
behaviors. Fig.4e and Fig.4f show that the new universal 
regime can be quite large. Finally we have also calculated 
the conductance fluctuation for 2D systems and similar 
behaviors were found (see Fig. 2d) In particular, the 
values of new UCF are found to be UCF(/3 = 2) = 0.68± 
0.02 and UCF(/3 = 4) = 0.66 ± 0.01. 

Now we provide further evidence of the universal be- 
havior of new UCF. In 2D systems there can be a metal- 
insulator transition (MIT) for P = 2, 4 but not for p = 1. 
This somehow coincides with our findings that there is 
a new UCF regime for /3 = 2,4 but not for P = 1. To 
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FIG. 5: (Color online) Localization length, average conduc- 
tance, conductance fluctuation of 2D symplectic systems for 
different L with E = 2.3 and t so = 0.4. (a). vs W. (b). 

rms(G) vs (G) for different L = 60, 70, 120. (c). rms(G) 
and (G) vs W. 



explore this correspondence further, we have calculated 
the localization length of the 2D system using the trans- 
fer matrix approach^ Here we calculate the localization 
length £ for quasi-lD systems with fixed length 1000000 
and different widths L\ (see Fig. 5a). Fig. 5a depicts the 
localization length vs disorder strength for (3 = 4. 

The intersection of different curves gives an estimate of 
the critical disorder strength of MIT of the 2D system. 
Fig. 5a shows that for an infinite 2D system, there is a 
MIT around W c — 6.1. For a mesoscopic 2D system, 
the critical region becomes a crossover region around the 
same W c and it is in this region where the new UCF is 
found. To see how our new UCF evolves with increas- 
ing of system size L we have calculated rms{G) vs W 
for finite 2D systems with different sizes L — 50 + lOn 
where n = 1, 2, 7. As shown in Fig. 5c, for a fixed W 
that is beyond the crossover regime, e.g., W = 7, the 
fluctuation decreases as L increases so that rms(G) — > 
at L —>■ oo. Importantly, the second plateau value (the 



new UCF) does not change with L. In addition, both 
(G) and rms{G) converge at W c — 6.1 for different L. 
This means that when L goes to infinity we should have 
rms(G) — cp in the vicinity of critical region where cp 
is the new UCF. This again suggests that the new UCF 
is driven by MIT and is an universal quantity. Finally 
rms(G) vs (G) is plotted in Fig. 5b for different L which 
shows the universal behavior that is also independent of 
L. Similar behaviors were also observed for 2D systems 
with /3 = 2. 

We have further calculated the third and fourth 
moments (71 and 72) of conductance vs the average 
conductance^ For QD and 2D systems similar univer- 
sal features found in rms(G) were also found for 71 and 
72 with different (3. 18 It is interesting that in the local- 
ized regime with (G) < 0.3 the conductance distribution 
seems to be superuniversal that is independent of dimen- 
sionality and symmetry^ 

In summary, we have carried out extensive simulations 
on conductance fluctuations of quasi-lD, QD, and 2D 
mesoscopic systems for orthogonal, unitary, and sym- 
plectic ensembles. Our results show that in addition to 
the usual UCF in the diffusive regime there exists a new 
UCF in the crossover regime between metallic and in- 
sulating regimes for unitary and symplectic ensembles 
but not for the orthogonal ensemble. We found that the 
conductance fluctuation rms(G) and higher order mo- 
ments vs (G) are universal functions from diffusive to lo- 
calized regimes which depend only on symmetry index (3 
for quasi-lD, QD, and 2D mesoscopic systems. In quasi- 
lD systems this universal function agrees with the result 
from DMPK equation. Our analysis suggests that this 
new UCF is driven by MIT in 2D systems. 
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In fig.2 of EPAPS, the rms(G) vs (G) is depicted in fig.2a. 
The conductance distribution for (G) = 0.2 shows univer- 



sal behavior independent of dimensionality and symmetry 
(fig. 2b) while for (G) = 0.6 deviation from universal be- 
havior is clearly seen (fig. 2c). 
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